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Abstract. Nuclear energy released by splitting Uranium and Thorium isotopes into 
two, three, four, up to eight fragments with nearly equal size are studied. We found that 
the energy released from equally splitting the 235 - 238 U and 230,232 Th nuclei into to three 
fragments is largest. The statistical multifragmentation model is employed to calculate 
the probability of different breakup channels for the excited nuclei. Weighing the the 
probability distributions of fragments multiplicity at different excitation energies for 
the 238 U nucleus, we found that an excitation energy between 1.2 and 2 MeV/u is 
optimal for the 235 U, 238 U, 230 Th and 232 Th nuclei to release nuclear energy of about 
0.7-0.75 MeV/u. 


| Email: imqmd@qq.com 
§ Corresponding author: fszhang@bnu.edu.cn 


Nuclear energy release from fragmentation 


2 


Nuclear energy is one of the most efficient sources of energy based on the Einstein’s 
mass-energy equivalence formula (A E = A me 2 ). It originates in the strong force holding 
the protons and neutrons together. The release of nuclear energy happens when nucleons 
rearrange themselves to form more stable nuclei. In such nuclear reaction processes, a 
large amount of energy is released in the form of emitted kinetic energy of particles or 
fragments and electromagnetic radiation. 

There are three ways to release nuclear energy, known as radioactive decay m, 
fusion [5TfT0] and fission [TT1 - 1I5] . The radioactive decay is the spontaneous process which 
occurs in radioactive materials by which the nuclei of unstable parent nuclei, gradually 
break up and are transformed into more stable isotopes or into nuclei of a different 
type. The daughter nuclei, consequently lose energy by emitting radiation in the form 
of particles (include a, f3, photons, neutrons particle etc.) and/or electromagnetic 
waves. The continuous fusion reaction happens in the core of stars including our 
sun. However in order to harness it to produce power, the fusion process must be 
controlled. Nuclear energy from fusion is still far from commercially viable even though 
controllable fusion reaction has been studied for many years through the means of 
magnetic confinement mm and inertial confinement PEH22]. The fission technology 
has been widely used to generate power by neutron-induced chain reactions. This 
fission process occurs when a heavy nucleus such as 235 U and 239 Pu absorbs a thermal 
neutron and the excited compound nucleus is excited to be near its fission barrier. 
At low excitation energy, quantum effects in the fission process are very significant 
to the dynamical evolution of the system. Nuclear structure effects, such as the shell 
corrections, directly influence the fission barrier height and the fission path. Most fission 
modes induced by thermal neutron are binary fissions in which one daughter nucleus has 
a mass of about 90 to 100 u and the remaining nucleus 130 to 140 u. Occasionally ternary 
fission and quaternary fission occur with relatively probability of 1/300 and 1/3000 per 
binary fission, respectively [231124]. For each binary event of uranium nucleus, about 0.1 
percent of the mass of the parent nucleus appears as the fission energy which is about 0.85 
MeV per nucleon. For the nucleus with high excitation energies, the multifragmentation 
characterized by the multiple production of intermediate mass fragments is the typical 
decay mechanisms which is quite different from that at low excitation energies. It is not 
yet clear how much energy will be released in the multifragmentation process. 

To describe the breakup process of excited nucleus, statistical models have been 
proposed, such as the statistical multifragmentation model (SMM), the statistical 
evaporation model (HIVAP) [251126] and the statistical model (GEMINI) [27]l28j. among 
others [29H3T] . The SMM used to describe the breakup of excited nuclei in this work has 
been successfully applied to describe intermediate energy heavy-ion collisions [32] and 
spallation reactions [33]. The basic assumption is that, after the most violent stages of 
the reaction, the source at thermal equilibrium with mass number Aq, charge number 
Z 0l and excitation energy E* will expand and desintegrate almost simultaneously into 
a number of small and intermediate mass fragments. For a given excitation energy E*, 
especially for highly excited nuclei, a large number of breakup channels are open. The 
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different breakup channels are generated by Monte Carlo simulations. The probability of 
each breakup channel is normally formulated in terms of a statistical partition function 
Wi which is related to the total energy and the Helmholtz free energy through the 
standard thermodynamic expression. The relative probability for a breakup channel 
can be written as P = IF,:/ ^ IF,;. And the average value of a physical observable O 
should be calculated through < O >= OfWi/ J^IR. For more detailed description 
of the SMM see Refs. [32TT36] . 

In the aftermath of the nuclear reaction, if the total mass of the final system is 
lower than the initial mass, the loss in mass which is called mass defect is converted 
into a large amount of energy given by (A E = A me 2 ). The difference of total binding 
energy before and after the reactions is defined as: 

A E= Y, E b( A v Z j)~ (1) 

j (final) {(initial) 

where Eb(Ai, Zf) denotes the binding energy of i-th nucleus. 

The specific binding energy (binding energy per nucleon) reflects the interactions 
of the nucleons inside the nucleus, especially the short range nuclear force and the long 
range Coulomb forces. Fig. 1 shows the experimental specific binding energy for 2438 
nuclei from the mass table AME2012 published in Refs. [37]. 56 Fe is the strongest 
bound nucleus and the nuclei with mass number around 56 have the largest specific 
binding energy values which is approximately 8.79 MeV/u. In general, for lighter nuclei 
with mass number less than 56, the nuclear force which is short range and attractive, 
has not reached saturation. Except for very light nuclei with increasing mass number, 
the average number of nucleons around one nucleon increases gradually, therefore, the 
specific binding energy also increases gradually. However, the Coulomb force, which is 
proportional to Z 2 , is long range and repulsive and become more important with heavier 
nuclei. The specific binding energy for heavy nuclei decreases gradually with increasing 
of mass number. 

In massive stars, the nuclear energy in a contracting core can be released through 
the neon, oxygen and silicon burning processes [38] • Complete exhaustion of the 
nuclear fuel in the stellar nucleosynthesis process happens when final production of 
the stable isotope 56 Fe is formed. In a typical fission process, such as the reaction 
of n+ 235 U—> 141 Ba+ 92 Kr+3n, the total energy released is about 0.85 MeV/u which 
including 0.70 MeV/u and 0.15 MeV/u corresponds to the energy released in the fission 
process and decay process of prefragments 141 Ba and 92 Kr, respectively. Since the mass 
number of 141 Ba and 92 Kr are much larger than that 56 Fe, the energy released in this 
fission process is incomplete. More nuclear energy can be released if such heavy nucleus 
can split into three or more fragments. Here we demonstrate the nuclear energy released 
for the three-fragments case. For simplicity, the decay process of unstable nuclei is not 
taken into account. We split a source with mass number A into three fragments A±, A 2 
and A 3 , and 


A — A\ + A 2 + A 3 , 


( 2 ) 
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where, A, = 0,1, 2, 3 • • • and each data point A is plotted in a Dalitz-type plot. For each 
set of fragments, the nuclear energy released by the source is calculated by Eq. (1) with 
mass table AME2012. Fig. 2(a) is a triangular Dalitz-type plot for 238 U. In a triangular 
Dalitz-type plot one can represent each event with fragments’ mass numbers A 1; A 2 , 
and A 3 by a point placed at distances from three sides of the equilateral triangle. A 
point located at the corner, on the side and inside the triangle, corresponds to an event 
with one, two, and three fragments, respectively. In the central region of Dalitz-type 
plot, the energy emitted by 238 U is 0.95 MeV/u as compared to 0.7~0.8 MeV/u released 
from splitting 238 U into two equal size fragments. If a heavy source such as 238 U splits 
equally into three fragments, the mass number of each fragments is about 80. To bring 
the fragments closer to 56 Fe, it is necessary to split it into four fragments. 

Similar to the ternary fragment case, we represent four fragment events with 
fragments’ mass numbers A 4 , A 2 , A 3 and A 4 by a point placed at distances from four 
surface of the tetrahedron. A 4 = 0 (Fig. 2(a)) denotes one of the four surface of the 
tetrahedron. Fig. 2(b) shows a slice of tetrahedron’s Dalitz plot calculated with the 
mass table AME2012 with A 4 = 60 for 238 U. From Fig. 2(b), one can see that 0.80 
MeV/u is released if the 238 U nucleus is splitted into four equal fragments. This is 
lower than 0.95 MeV/u obtained by splitting the nucleus into three equal fragments, 
despite the mass number of each of the four fragments be about 60 which is closer to 56. 
The fragments splitted from source 238 U have larger neutron-proton ratio, the specific 
binding energy of these fragments is very small as the symmetry energy term plays an 
important role in nuclear binding energy. 

Fig. 3 compares the average energy released per nucleon obtained with the mass 
table AME2012 for splitting 235 U, 238 U, 230 Th and 232 Th into two, three, four, five, six, 
seven, and eight fragments with nearly equal size. Nearly equal energy is released 
in the case of the two, four, and five fragments, which is about 0.8 MeV/u. The 
energy released from three fragments case is larger than in the other cases. Although 
splitting it into three fragments with nearly equal size is the most efficient way of 
releasing nuclear energy, the probability of ternary fission for the heavy nucleus at 
low excitation energies is very small. It is necessary to increase the excitation energy 
of the source to study the nuclear energy released in the multifragmentation process. 
However the multifragmentation process for highly excited nuclei is very complex, as 
a large number of decay channels are simultaneously open [39]. We employ statistical 
multifragmentation model (SMM) which is a classical statistical model to calculate the 
probability of different decay channels for the excited nucleus. 

Fig. 4 shows the probability as a function of the fragment multiplicity for 238 U 
sources with excitation energy E* =1, 2, 3 and 4 MeV/u. We select events whose 
multiplicity values lie between 2-10. It is seen that many breakup channel are 
simultaneously open for a given excitation energy. For E* = 1 MeV/u, the probability 
are approximately 0.6 for binary events and 0.4 for ternary events. As the excitation 
energy increases, decay channels with larger multiplicities are opened gradually and the 
position of the maximum probability is moves to higher multiplicity values. 
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In the framework of SMM, the average energy released is obtained by averaging this 
quantity in each fragment aion mode < A E >= ATyllq/ ]C W, n where AE, denotes 
the nuclear energy release in the i-th event and Wi is corresponding statistical weight. 
The average nuclear energy released per nucleon for the 235 U, 238 U, 230 Th, and 232 Th 
at different excitation energies are shown in Fig. 5. As is well known, for a heavy 
nuclei with low excitation energy, one observes the emission of light particles plus an 
evaporation residue in the deexcitation process. Fission occurs usually when the excited 
energy is larger than the fission barrier of the nuclei. Multifragmentation often happens 
in highly excited nuclei. From Fig. 5 one can see that at E* =1 MeV/u, the average 
nuclear energy release rises up to 0.7 MeV/u, as one varies the source nucleus, due to the 
opening of the binary breakup channel. Between 1.2 and 2 MeV/u, the probability of 
binary, ternary and quaternary breakup channel accounts for the largest proportion of 
the emitted fragments, and peaks around 0.7~0.75 MeV/u. However as the excitation 
energy surpasses 2 MeV/u, the breakup channels with higher multiplicities dominates 
the partition function. The average energy released per nucleon gradually decreases 
with the increasing of the excitation energy. 

In summary, the nuclear energy released has been calculated by splitting 235 U, 
238 U, 230 Th and 232 Th sources into two, three, four, five, six, seven, and eight fragments 
with nearly equal size. The calculation results demonstrate that fragmentation into 
three nearly equal fragments is the most efficient mechanism to release the nuclear 
energy. We found that 0.95 MeV per nucleon is released in this case. The statistical 
multifragmentation model has been employed to simulate the breakup process of 
excitation nucleus. Taking into account the probability distributions of fragments 
multiplicity at different excitation energies we ford that an excitation energy between 
1.2 and 2 MeV/u is optimal for the 235 U, 238 U, 230 Th and 232 Th nuclei to release nuclear 
energy. 

Acknowledgments 

We wish to thank J. Randrup for fruitful discussion. This work was supported by 
the National Natural Science Foundation of China under Grants No. 11025524 and 
No. 11161130520 and the National Basic Research Program of China under Grant No. 
2010CB832903. 

References 

[1] Behrens H and Biihring W 1971 Nucl. Phys. A 162 111. 

[2] Holstein B R, Shanahan W and Treiman S B 1972 Phys. Rev. C 5 1849. 

[3] Childress J D 1961 Phys. Rev. 123 1729. 

[4] Rose M E and Osborn R K 1954 Phys. Rev. 93 1326. 

[5] Jones S E 1986 Nature 321 127. 

[6] Pease R S 1958 Nature 182 1051. 

[7] Lawson J D 1957 Nature 180 780. 



Nuclear energy release from fragmentation 


6 


[8] Pampin R, Davis A, Izquierdo J, Leichtle D, Loughlin M J, Sanz J, Turner A, Villari R and Wilson 

P P H 2013 Fusion Eng. Des. 88 454. 

[9] Sawyer R F 2010 Phys. Rev. Lett. 104 191103. 

[10] Lawson J D 1957 Proc. Phys. Soc. B 70 6. 

[11] Meitner L and Frisch O R 1939 Nature 143 239. 

[12] Hahn O and Strassmann F 1939 Naturwiss. 27 11. 

[13] Kruger P and Sugarman N 1955 Phys. Rev. 99 1459. 

[14] Lavrukhina A K and Krasavina L D 1957 Atomic Energy 2 27. 

[15] Diven B C, Terrell J and Hennnendinger A 1958 Phys. Rev. 109 144. 

[16] Strait E J 2015 Phys. Plasmas 22 021803. 

[17] Sun X Y, Wang F, Wang Y and Li S 2014 Fusion Engineering and Design 89 741. 

[18] Wan B N, Li J G, Guo H Y, Liang Y F, Xu G S and Gong X Z 2013 Nucl. Fusion 53 104006. 

[19] Jafari S, Nilkar M, Ghasemizad A and Mehdian H 2014 Phys. Plasmas 21 104503. 

[20] Bychkov V, Modestov M and Law C K 2015 Progress in Energy and Combustion Science 47 32. 

[21] Eder D C, Fisher A C, Koniges A E and Masters N D 2013 Nucl. Fusion 53 113037. 

[22] Li C K, Seguin F H, Frenje J A, Rosenberg M J, Rinderknecht H G, Zylstra A B, Petrasso R D, 

Amendt P A, Landen O L, Mackinnon A J, Town R P J, Wilks S C, Betti R, Meyerhofer D D, 
Soures J M, Hund J, Kilkenny J D and Nikroo A 2012 Phys. Rev. Lett. 108 025001. 

[23] Tsien S T, Ho Z W, Vigneron L and Chastel R 1947 Nature 159 773. 

[24] Titterton E W 1952 Nature 170 794. 

[25] Wang N, Zhao K, Scheid W and Wu X Z 2008 Phys. Rev. C 77 014603. 

[26] Vinodkumar A M, Loveland W, Sprunger P H, Prisbrey L, Trinczek M, Dombsky M, Machule P, 

Kolata J J and Roberts A 2009 Phys. Rev. C 80 054609. 

[27] Charity R J 2010 Phys. Rev. C 82 014610. 

[28] Charity R J, Sobotka L G, Cibor J, Hagel K, Murray M, Natowitz J B, Wada R, Masri Y El,Fabris 

D , Nebbia G, Viesti G, Cinausero M, Fioretto E, Prete G, Wagner A and Xu H 2001 Phys. 
Rev. C 63 024611. 

[29] Das C B, Gupta S Das, Lynch W G, Mekjian A Z and Tsang M B 2005 Phys. Rep. 406 1. 

[30] Moretto L G and Wozniak G J 1993 Annu. Rev. Nucl. Part. Sci. 43 379. 

[31] Bondorf J P, Botvina A S, Iljinov A S, Mishustin I N and Sneppen K 1995 Phys. Rep. 257 133. 

[32] D’Agostino M, Botvina A S, Bruno M, Bonasera A, Bondorf J P, Bougault R, Desesquelles P, 

Geraci E, Gulminelli F, Iori I, Neindre N L, Margagliotti G V, Mishustin I N, Moroni A, Pagano 
A and Vannini G 1999 Nucl. Phys. A 650 329. 

[33] Abdel-Waged K, Felemban N, Gaitanos T, Ferini G and Toro M Di 2010 Phys. Rev. C 81 014605. 

[34] Donangelo R, Sneppen K and Souza S R 2001 Comput. Phys. Commun. 140 405. 

[35] Tan W P, Souza S R, Charity R J, Donangelo R, Lynch W G and Tsang M B 2003 Phys. Rev. C 

68 034609. 

[36] Souza S R, Carlson B V, Donangelo R, Lynch W G and Tsang M B 2013 Phys. Rev. C 88 014607. 

[37] Wang M, Audi G, Wapstra A H, Kondev F G, MacCormick M, Xu X and Pfeiffer B 2012 Chinese 

Phys. C 36 1603. 

[38] Heger A, Woosley S E, Rauscher T, Hoffman R D and Boyes M M 2002 New Astronomy Reviews 

46 463. 

[39] Sneppen K 1987 Nucl. Phys. A 470 213. 



Nuclear energy release from fragmentation 


7 



0 _i_i_i_i_i_i_i_i_i_i_t 

0 50 100 150 200 250 


A 


300 


Figure 1 . Empirical specific binding energy as a function of mass number A for 2438 
nuclei from the new mass table AME2012 m- 
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Figure 2. (a) The triangular Dalitz-type plot for the nuclear energy released obtained 
with empirial values by splitting 238 U into three fragments. The mass numbers of the 
three fragments and the source are denoted by A i: (i = 1,2,3) and A, respectively, 
(b) The slice of tetrahedron’s Dalitz plot for 238 U with A 4 = 60. The A' denotes 
the distances between one point and one of three side on slice and satisfy relation 
A[ = 3\/2/ 4Ai (i = 1, 2, 3) which is obtained by geometric relationships. 
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Figure 3. Average energy released per nucleon for 235 U, 238 U, 230 Th and 232 Th split 
to 2, 3, 4, 5, 6, 7, and 8 fragments with nearly equal size. 



Figure 4. The normalized multiplicity distribution! as a function of fragment 
multiplicity obtained with SMM for E* =1, 2, 3 and 4 MeV/u. 
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Figure 5. Average energy release as a function of the excitation energy E* = 1,2,3 
and 4 MeV/u for the 2 32,238 Th and 235.238U sources . 






